Stability 



1. The Notion of Stability 

A system is called stable if its long-term behavior does not depend sig- 
nificantly on the initial conditions. 

It is an important result of mechanics that any system of masses con- 
nected by springs (damped or undamped) is a stable system. In network 
theory there is a similar result: any RLC-network gives a stable system. In 
these notes, we investigate for the simplest such systems why this is so. 

In terms of differential equations, the simplest spring-mass system or 
RLC-circuit is represented by an ODE of the form 

floy" + «iy' + «2y = ^(0' constants, t = time. (1) 

For the spring-mass system, y is the displacement from equilibriimi po- 
sition, and r{t) is the externally applied force. 

For the RLC-circuit, y represents the charge on the capacitor, and r{t) is 
the electromotive force £{t) applied to the circuit (or else y is the current 
and r(0 = 

By the theory of inhomogeneous equations, the general solution to (1) 
has the form 

y = Ciyi + C2y2 + yp, C\, Ci arbitrary constants, (2) 

where y^ is a particular solution to (1), and Ciyi + C2y2 is the complementary 
function, i.e., the general solution to the associated homogeneous equation 
(the one having r{t) =0). 

The initial conditions determine the exact values of ci and ca. So from 

(2), 

the system modeled for every choice ofci, C2, 

by (1) is stable ciy\ + C2y2 — >■ as i -> oo. 

Often one applies the term stable to the ODE (1) itself, as well as to the 
system it models. We shall do this here. 

If the ODE (1) is stable, the two parts of the solution (2) are named: 

yp = steady-state solution c\y\ + C2y2 = transient; (4) 
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the whole solution y{t) and the right side r{t) of (1) are described by the 
terms 

y{t) = response r{t) = input. 

From this point of view, the driving force is viewed as the input to the 
spring-mass system, and the resulting motion of the mass is thought of as 
the response of the system to the input. So what (2) and (4) are saying is that 
this response is the sum of two terms: a transient term, which depends on 
the initial conditions, but whose effects disappear over time; and a steady- 
state term, which represents more and more closely the response of the 
system as time goes to oo, no matter what the initial conditions are. 

2. Conditions for Stability: Second Order Equations 

We now ask under what circumstances the ODE (1) will be stable. In 
view of the definition, together with (2) and (3), we see that stability con- 
cerns just the behavior of the solutions to the associated homogeneous 
equation 

aoy" + flij/' + fl2i/ = ; (5) 

the forcing term r{t) plays no role in deciding whether or not (1) is stable. 

There are three cases to be considered in studying the stability of (5); 
they are summarized in the table below, and based on the roots of the char- 
acteristic equation 

aov'^ + a\r + = Q ■ (6) 

roots solution to ODE condition for stability 

r\ n cie''!* -f C2e''2* ri < 0, r2 < 

ri = r2 e''i*(ci -I-C2O ri < 

a ± ih e"* {c\ cos ht -\- c^?>m.ht) a < 

The first two columns of the table should be familiar, from your work in 
solving the linear second-order equation (5) with constant coefficients. Let 
us consider the third column, therefore. In each case, we want to show that 
if the condition given in the third column holds, then the criterion (3) for 
stability will be satisfied. 

Consider the first case. If ri < and r2 < 0, then it is immediate that 
the solution given tends to as i 00. 

On the other hand, if say r\ > 0, then the solution e*"!* tends to 00 (or to 1 
if fi =0). This shows the ODE (5) is not stable, since not all solutions tend 

to as f — )■ 00. 
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In the second case, the reasoning is the same, except that here we are 
using the limit 

lim te'^ = ^ r <0 

For the third case, the relevant limits are (assuming b 7^ for the second 
limit): 

lim e"'^ cosbt = <^ a < 0, lim e"^ smbt = <^ a <0 . 

The three cases can be summarized conveniently by one statement: 
Stability criterion for second-order ODE's — root form 

// , / , • all roots of Unr^ + air + ao = 

any + aiy -\-a2y = rlt) is stable ^ , , . , , 

nave negative real part. 

(7) 

Alternatively one can phrase the criterion in terms of the coefficients of the 
ODE; this is convenient, since it doesn't require you to calculate the roots 
of the characteristic equation. 

Stability criterion for second order ODE's — coefficient form. Assume 
flo > 0. 

aoy" + fliy' + a2y = r{t) is stable <^=^ flo/fli/«2 > . (8) 
The proof is left as an exercise; it is based on the quadratic formula. 
3. Stability of Higher Order ODE's 

The stability criterion in the root form (7) also applies to higher-order 
ODE's with constant coefficients: 

(aoD" + aiD"-^ + ... + a„_iD + fl„)y = f{t). (9) 

These model more complicated spring-mass systems and multi-loop RLC 
circuits. The characteristic equation of the associated homogeneous equa- 
tion is 

aor'' + air""'^ + . . . + an-ir + an = 0. (10) 

The real and complex roots of the characteristic equation give rise to solu- 
tions to the associated homogeneous equation just as they do for second 
order equations. (For a /c-fold repeated root, one gets additional solutions 
by multiplying by 1, t, t^,. . . i'^~^.) 
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The reasoning which led to the above stability criterion for second- 
order equations applies to higher-order equations just as well. The end 
result is the same: 



Stability criterion for higher-order ODE's — root form 

ODE (9) is stable <^=^ all roots of (W) have negative real parts; 



(11) 



that is, all the real roots are negative, and all the complex roots have nega- 
tive real part. 

There is a stability criterion for higher-order ODE's which uses just the 
coefficients of the equation, but it is not so simple as the one (8) for second- 
order equations. We will not use this in the course, but it is worth seeing. 
The key point is that the stability of the system can be found without find- 
ing the roots of a higher order poljmomial. 

Without loss of generality, we may assume that aq > 0. Then it is not 
hard to prove that 



ODE (9) is stable 



ao,...,a„ > . 



(12) 



The converse is not true. For an implication <^=, the coefficients must sat- 
isfy a more complicated set of inequalities, which we give without proof, 
known as the 



Routh-Hurwitz conditions for stability 

stable ^ 



Assume ao > 0; ODE (9) is 



in the determinant below, all of the n principal minors (i.e., the subdeterminants in 
the upper left corner having sizes respectively 1,2, ... ,n) are > when evaluated. 



(13) 



In the determinant, we define a^^ = Q if k > n; thus for example, the last row 
always has just one non-zero entry, an. 
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